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Abstract 

An  orbiting  detector  of  Infiared  (IR)  en»gy 
may  be  used  to  detect  the  rocket  plumes  generated 
by  ballistic  missiles  during  the  powoed  segment 
of  their  trajectory.  By  measuring  angular 
directions  of  the  detections  over  several 
observations,  the  trajectory  propoties,  launch 
location,  and  impact  area  may  be  estimated  using 
a  nonlinear  least-squares  iteration  procedure. 
Observations  from  two  or  more  sensors  may  be 
combined  to  form  stereoscopic  lines-of-sight 
(LOS)  increasing  the  accuracy  of  the  algorithm. 
Presented  in  this  paper  is  a  computer  model  of  an 
estimation  algorithm  which  determines  what 
parameter,  or  combination  of  parameters,  will 
have  a  significant  effect  on  the  error  of  the 
tactical  parameter  estimation.  This  model,  which 
is  coded  in  MATLAB™,  generates  observation 
data  and  then,  using  the  data,  produces  an 
estimate  of  the  tactical  parameters,  i.e.,  the  time, 
position,  and  heading  at  launch  and  burnout,  and 
an  impact  time  and  position.  The  effects  of  time 
eiTOTS,  LOS  measurement  errors,  and  satellite 
position  errors  on  the  estimation  accuracy, 
determined  using  analytical  and  Monte  Carlo 
simulation  techniques,  are  presented. 

Introduction 

The  TRW-built  Defense  Support  Program 
(DSP)  satellites  have  been  the  spacebomib 
segment  of  NORAD’s  Tactical  Warning  and 
Attack  Assessment  system  since  the  early 
1970’s.  Using  infrared  detectors  that  sense  the 
heat  from  missile  plumes  against  the  Earth 
background,  these  orbiting  sentries  detect 
ballistic  missile  launches.  The  DSP  system 
provides  near  real  time  detection  information  in 
support  of  DOD’s  tactical  warning  and  attack 
assessment  mission,  and  is  supported  by  a 
network  of  fixed  and  mobile  ground  stations  that 


process  and  disseminate  information  to  military 
commanders  worldwide.  The  Cold  War  mission 
of  the  DSP  system  was  to  detect  massive 
intercontinental  ballistic  missile  (ICBM)  attacks. 
The  United  States’  response  to  such  an  attack 
only  required  timely  and  ambiguous  warning  and 
missile  flight  times  were  much  longer  than  the 
time  required  to  launch  a  retaliatory  attack. 
Precise  radar  tracks  could  be  established  with 
enough  time  to  mitigate  effects  as  much  as 
technology  allowed. 

The  current  combat  environment  demands 
much  more  from  launch  detection  satellites.  The 
present  threat  is  from  tactical  missiles  (TBMs), 
which  exhibit  much  cooler  and  shorter  thrust 
times,  and  possibly  more  dqnessed  trajectories 
than  those  exhibited  by  ICBMs.  TBMs  can  be 
launched  from  almost  anywhere  within  a  large 
geographical  area  of  interest,  with  lofted  or 
dejxes^  trajectories.  To  attack  the  launchor 
and/or  employ  anti-ballistic  missile  (ABM) 
weapons  or  alert  potential  victims  within  the 
impact  zone  the  tactical  parametras  must  be 
estimated  faster  and  more  accurately  than  was 
required  for  the  massive  ICBM  launch.  For 
budgetary  reasons  the  U.S.  is  forced  to  use  the 
existing  DSP  system  to  countw  the  TBM  threat 
into  the  beginning  of  the  next  century*  until  the 
Space  Based  Infrared  System  (SBIRS)  is 
operational.  More  rapid  extraction  of  more 
precise  information  from  these  existing, 
technology  limited  satellites  must  be 
accomplished  in  the  interim.  Anti-ballistic 
missile  (ABM)  systems,  such  as  Patriot  or 
Aegis,  may  be  employed  as  ABM  umbrellas  in  a 
tactical  area  of  interest,  provided  they  receive 
precise  and  timely  cueing  from  DSP. 

This  raises  the  question:  “How  accurate  is 
the  information  provided  by  DSP?”  To  answCT 
this  question  an  engineering  error  analysis  was 
performed  to  determine; 

•  What  are  the  errors  present  in  the  detection 

system? 


RMfEK 


QVALITV 


1 


•  Which  errors  have  the  greatest  effect  upon 

the  accuracy  of  DSP  output  information? 

•  What  are  the  effects  of  the  errors,  both 

individually  and  collectively,  on  the 

estimated  launch  and  impact  points? 

These  questions  are  addressed  in  this  paper.  The 
algorithm  used  by  the  tactical  warning  system 
(TALON  SHIELD/ALERT)  to  determine  the 
missile  trajectcxy  and  predict  the  launch  and 
impact  points  is  modeled,  the  various  errors  are 
int^uc^  and  their  effect  is  detomined  both 
analytically  (statistically)  and  numerically  by 
simulation.  The  results  presented  in  this  p^)ct 
are  taken  from  the  Masters  thesis  of  the  first 
authoP. 

This  paper  differs  from  the  earlier  work  by 
Danis^  who  considered  the  estimation  of  the 
launch  point  and  trajectory  from  just  two 
observations.  In  addition,  in  this  paper  the  errors 
in  impact  time  and  point  are  analyzed  and  the 
approach  is  different. 

DSP 

The  DSF*  system  consists  of  one  or  more 
satellites  in  geosynchronous  orbit  and  one  or 
more  ground  receiving  stations.  The  DSP 
satellite  is  10  meters  long,  7  meters  in  diameter, 
and  weighs  over  2300  kilograms  (Fig.  1).  Each 
satellite  is  spun  at  six  rpm  along  its  longitudinal 
axis  with  the  telescope  pointing  toward  the  Earth 
so  that  the  sensor  scans  the  Earth.  A  counter¬ 
rotating  wheel  keeps  the  system  in  a  nominal 
zero  momentum  state.  This  configuration  is 
sometimes  called  a  “yaw  spinner”.  The  IR 
telescope  is  tilted  from  the  spin  axis,  so  that  the 
photo-electric  cell  (PEC)  array  covers  the  radius 
of  the  Earth.  As  the  satellite  rotates,  the  entire 
surface  of  the  Earth  within  the  field  of  view 
(FOV)  is  scanned  by  the  IR  detector  as  shown  in 
Figs.  2  and  3. 

Detection  of  IR  sources  is  accomplished 
with  the  telescope  and  PEC  array  portions  of  the 
sensor.  The  PEC-array  of  the  IR  detector  is 
mounted  with  the  nadir  end  at  the  center  of  the 
rotation  of  the  telescope  (Fig.  2).  The  array 
contains  over  6000  detector  cells  that  are 
sensitive  to  energy  in  the  infrared  wavelengths. 
As  the  PEC  array  scans  the  FOV,  a  cell  passing 
across  an  IR  source  will  generate  a  voltage  with 
an  amplitude  proportional  to  the  signal  intensity. 
This  voltage  signal  is  termed  an  IR  return,  and  is 


transmitted  to  ground  processing  stations  after 
amplification  and  background  filtering.  Line  of 
sight  of  the  IR  source  relative  to  the  satellite  is 
detomined  from  the  angle  of  rotation  of  the 
telescope  at  the  time  of  detection  and  the  cell 
illuminated. 

Error  Analysis  Algorithm 

The  quality  of  the  trajectory  estimation 
process  is  of  paramount  importance  to  Tactical 
Ballistic  Missile  Defense  (TBMD).  Real-time 
knowledge  of  the  launch  position  allow  targeting 
of  the  launcher.  Cueing  for  ABM  systems,  such 
as  Patriot  and  Aegis,  requires  timely  and  accurate 
trajectory  information,  which  can  be  propagated 
from  knowledge  of  the  state  vector  at  burnout. 
Impact  time  and  position  is  extrapolated  from  the 
state  vector  at  burnout,  and  may  be  used  for 
warning  personnel  within  the  target  area. 
Understanding  the  algorithms  and  equations 
employed  in  the  estimation  process  is  necessary 
to  assess  the  quality  of  the  estimated  parameters. 

The  detection  of  a  TBM  launch  is  the 
starting  point  for  any  ballistic  missile  defense. 
DSP  does  this  by  detecting  IR  radiation  emitted 
by  the  exhaust  plume  of  a  launching  missile. 
With  detections  by  two  or  more  spacecraft 
(stereo)  ,  triangulation  of  lines  of  sight  can  be 
used  to  more  accurately  estimate  the  boost-phase 
trajectory,  which  is  then  used  to  calculate  launch 
position,  state  vector  (position  and  velocity)  at 
engine  burnout,  and  impact  position. 

The  tactical  parameter  estimation  process  is 
composed  of  several  tasks:  Initial  estimate  of  the 
tactical  parameters,  nonlinear  least-squares 
estimation  (refinement)  of  the  tactical  parameters, 
burn-out  time  estimation,  state  vector  generation, 
and  impact  point  calculation.  The  tactical 
parameters  are: 

Tq  =  Time  of  launch 
L  =  Loft 

00  =  Launch  point  geodetic  latitude 
Ao  =  Launch  point  geodetic  longitude 
Hq  =  Launch  point  height  above  WGS  -  84 

ellipsoid 

ao  =  Flight  trajectory  azimuth  (true  heading) 
Observational  Data 

From  the  observational  data  the  LOS  relative 
to  the  satellite  is  determined.  With  the  attitude 
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Table  1  Example  Observational  Data 


Index 

S/C 

Intens 

Azimuth 

Elevation 

G.H.A. 

Dec. 

Radius 

ID 

(rad) 

(rad) 

(rad) 

(rad) 

JM _ 

k 

S/C 

I 

lOHHHI 

11 

^ 

5 

R 

1 

129.36 

1 

14.0 

4.061854 

0.115847 

■iIlZUVJH 

0 

42164.17 

2 

2 

23.0 

2.427020 

0.094741 

1.221730 

0 

42164.17 

3 

135.44 

3 

32.0 

2.062181 

0.142310 

1.832595 

0 

42164.17 

4 

139.36 

1 

29.0 

4.062497 

0.115971 

0.174532 

0 

42164.17 

5 

140.30 

2 

29.0 

2.427264 

0.094753 

1.221730 

0 

42164.17 

6 

145.44 

3 

40.0 

2.061969 

0.142405 

1.832595 

0 

42164.17 

DBi 

149.36 

1 

49.0 

4.063552 

0.116147 

0.174532 

0 

42164.17 

8 

150.30 

2 

60.0 

2.427660 

0.094746 

1.221730 

0 

42164.17 

9 

155.44 

3 

65.0 

2.061752 

0.142493 

1.832595 

0 

42164.17 

data  which  is  included  in  the  telemetry  and  the 
satellite  position  as  detamined  by  the  Air  Force 
Satellite  Control  Network  (AFSCN)  the  LOS 
from  the  satellite  to  the  IR  event  in  Earth 
centaed  coordinates  can  be  determined.  Table  1 
shows  a  typical  set  of  observations  of  a  single 
TBM  launch. 

The  index,  k,  runs  from  1  to  n,  the  total 
number  of  observations,  and  is  used  as  a 
subscript  for  the  remaining  symbols  to  denote  a 
particular  observation.  Time,  T^,  is  the  time  of 
observation  in  seconds  measured  from  midnight 
of  the  day  of  the  observation.  Spacecraft 
Identification,  S/C,  identifies  which  satellite  is 
making  the  observation.  In  this  example  the 
event  was  detected  by  three  satellites.  The 
intensity,  1^,  is  the  radiant  intensity  of  the  IR 
return,  and  is  used  to  select  the  type  of  TBM 
being  detected.  Azimuth  angle,  is  the 
azimuth  LOS  of  the  IR  return  measured 
clockwise  from  true  south,  and  has  a  range  of 
(0,2ji)  radians  (Fig.  4).  Elevation  angle,  Tj^,  is 
the  elevation  of  the  return  measured  from  nadir 
and  has  values  between  (0,  0.175)  radians  or 
(0,10)  degrees.  The  satellite  position  is  given  in 
spherical  coordinates  (Fig.  5).  Greenwich  hour 
angle,  gha^,  is  the  angle  between  the  satellite’s 
nadir  point  and  the  prime  meridian  measured  east. 
Declination  angle,  5^,  is  the  angle  above  or 
below  the  equator  measured  positive  north. 
Radius,  R^,  is  the  distance  from  the  satellite  to 
the  Earth’s  center,  measured  in  kilometers. 


TBM  Profile 

A  TBM  profile  is  a  description  of  the 
nominal  powoed  flight  trajectory  of  a  given 
TBM.  A  profile  consists  of  IR  intensity  and  the 
nominal  (maximum  range)  vertical  and  horizontal 
ranges  from  the  launch  point  as  a  function  of 
time.  The  detected  radiant  intensity  is  ccxnpared 
to  TBM  profiles  in  a  data  base,  and  the  best 
match  is  selected  as  the  type  of  TBM  being 
observed.  This  selection  process  is  complex  and 
will  be  assumed  to  have  been  done  correctly  in 
this  study.  The  downrange  and  altitude  profiles 
are  rqwesented  by  fourth  order  polynomials  in 
time  of  the  form: 

dp=aQ  +  ait+ ,  downrange 
hp  =  bo  +  bit  +  b2t^ +  +  b4t*,  altitude  (1) 

where  t  is  the  time  of  flight  from  launch  and  the 
3;  and  b;  are  coefficients.  These  coefficients  ate 
nol  determined  by  the  algorithm,  but  are  assumed 
to  be  known  fiom  the  typing  process,  i.e., 
perfect  a  priori  knowledge  of  the  particular 
observed  TBMs  nominal  trajectory.  Since  real- 
life  TBMs  do  not  fly  the  trajectory  exactly  using 
the  profile  introduces  an  error.  This  error  is 
ignored  in  this  analysis  since  the  purpose  ho'e  is 
to  determine  the  effects  of  other  error  sources. 

A  loft  parameter,  L,  is  used  to  account  fw 
trajectories  above  or  below  the  nominal  profile. 
The  actual  range  and  altitude  are  modeled  by 

d  =  il-l.5L)dp 

h  =  (l  +  L)hp 
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The  loft  varies  over  ±0.25,  L>0  represents  a 
lofted  trajectory  and  L<0  represents  a  dqxessed 
trajectory. 

Nonlinear  Least  Squares  Estimation 

For  the  estimation  process  it  is  convenient 
to  express  the  “polar”  focal  plane  observations 
(/J.jj)  into  focal  plane  “Cartesian”  coordinates 
(k.v)  so  that  the  coordinates  have  similar 
behavior  with  respect  to  errors  and  noise.  The 
transformation  is 


H  =  -tanTjsinj3 
v  =  -tanT;cos)3 


(3) 


!2=(yo-ycrw(yo-yc) 

Q  =  {Sy-ASxfw{Sy-ASx)  (6) 

5y  =  (yo-yc) 


where  W  is  a  weight  matrix  and  is  usually  taken 
as  a  diagonal  matrix  where  the  diagonal  terms  are 
in  the  inverse  of  the  variances  of  the 
observations.  In  this  analysis  we  have  assumed 
W  is  the  identity  matrix.  We  now  want  to  find 
the  change  in  the  value  of  the  tactical  parameters, 
5x,  that  will  minimize  Q.  This  is 
accomplished  by  taking  the  partial  derivative  and 
setting  it  equal  to  zero. 


For  the  nonlinear  least  squares  estimation  let  y 
be  the  vector  of  observations,  and  let  X  be  the 
tactical  parameter  vector,  that  is 


y^={yi.y2*— *y„).yt  =K.vt) 


-^  =  2A'^W(5y-A&)  =  0  (7) 

dox 

If  there  are  n  observations  then  this  is  a  set  of  2n 
equations  with  six  unknowns.  The  solution  is 

8x  =  (  A'*'Wa)'^  A'^Wdy  (8) 


Note  that  vectors  are  represented  by  bold  fece 
letters,  yo  will  represent  the  actual  observations 
and  yc  the  calculated  observations  based  on 
specific  values  of  the  tactical  parameters.  The 
nonlinear  least  squares  proc^ds  as  follows: 
Assume  some  initial  values  of  the  tactical 
parameters,  Xq,  (how  to  obtain  these  will  be 
discussed  later)  expand  the  equations  for  the 
observations  in  a  Taylor  series  about  these  initial 
values,  and  retain  the  first  term  of  the  series. 


yo(x.0=y(xc»0+-^| 


(x-xc) 


yo(x,0=yc+A5x  (5) 

A  =  ^ 
ax__ 


A  is  the  partial  dmvative  of  the  observations 
with  respect  to  the  tactical  parameters.  In  the 
least  squares  approach  we  find  the  values  of  the 
tactical  parameters,  X,  that  minimize  the 
diff^ence  between  the  calculated  observations, 
yc,  and  the  actual  observations,  yo-  That  is, 
we  minimize 


The  process  is  shown  in  Figure  6.  With  an 
initial  estimate  of  X,  calculate  the  observations 
y^  and  the  matrix  of  partials  A,  then  with 
Sy  =  yo-yc  calculate  Sx.  (The  initial 
estimate  and  A  are  addressed  in  the  following 
two  sections.)  With  (xc  )j  =  (xc  )j_j  +  Sx 

continue  the  process  until  the  changes  in  Sx 
become  smaller  than  some  tolerance.  The 
resulting  solution  is  the  best  estimate  in  the  least 
squares  sense. 

If  W  is  the  inverse  of  the  covariance  of  the 
measurement  errors  then  it  can  be  shown  that  the 
matrix 

C  =  (a’'Wa)"‘  (9) 

is  the  covariance  matrix  of  the  tactical 
parameters. 

Initial  Estimate 

The  initial  guess  of  the  launch  point  is 
obtained  by  projecting  the  LOS  for  the  first 
observation  for  each  of  the  satellites  detecting  the 
launch  onto  the  surface  of  the  Earth  which  yields 
a  latitude  and  longitude.  These  values  are  then 
avCTaged  to  obtain  the  initial  guess.  For 
example,  if  the  first  three  observations  are  from 


4 


(16) 


diff^nt  satellites  and  0,-  and  A,-, i  =  1,2,3,  are 
the  latitudes  and  longitudes  of  the  projections  of 
the  observations  onto  the  surface  of  the  Earth, 
then  the  initial  estimates  are 


00°^  =(01  +02  +  03)/ 3 
4‘’^=(Ai  +  ;.2  +  A3)/3 


(10) 


The  initial  estimates  for  LandAo  ate 
L  =  0,  Ao  =  0 .  The  time  of  launch  is  taken  as  20 
seconds  before  the  first  observation,  i.e., 

ro  =  ri-20  (11) 


A.  —  A4A3A2AJ 

A4  is  the  transformation  of  the  obsCTvations 
from  the  focal  plane  into  the  satellite  Cartesian 
frame,  A3  then  transforms  them  into  the  Earth 
centered  Cartesian  frame  and  A2  transforms  them 
into  Earth  caiteed  spherical  coordinates.  Aj  is 
the  matrix  of  partial  derivatives  of  the  ^herical 
coordinates  with  respect  to  the  tactical 
parameters.  This  decomposes  A  into  four 
matrices  which  are  relatively  easy  to  daive 
analytically.  They  are  given  in  the  Appendix. 
The  reado"  is  refered  to  Ref.  2  for  their 
derivation. 


To  obtain  the  heading  the  last  observations  from 
each  of  the  satellites  are  avoaged  in  the  same 
manner  as  the  first  to  obtain 

4>last  =  {0«-2  +  0n-l  +  ^ 

^last  =  {^n-2  +  ^n-l  +  ^n)  /  3 

The  heading  is  obtained  from  the  great  circle  arc 
between  these  two  points. 


Burnout  Time  Estimation 

The  estimation  of  burnout  time  is  based 
on  the  TBM  profile  maximum  bum  time  (tmix)* 
the  last  observation  time  {T^),  and  the  next 
potential  observation  time  (T^xc)  I'^d  it 
occurred.  The  maximum  bum  time  according  to 
the  profile  is 

7’m.x  =  ro  +  W  (17) 


ar=f-tan-‘ 


[0)1 
“  J 

|COS0^°^^ 

1 

^  J 

(13) 


Two  cases  can  occun 

a)  If  r,„,>r„^then 


Matrix  of  Partials 

The  matrix  of  partials  is  best  developed  by 
decomposing  it  into  a  set  of  matrix 
multiplications.  Referring  to  Figure  5  the 
(X,Y2)  coordinate  is  fixed  to  the  Earth  and  has 
its  origin  at  the  Earth’s  center.  The  (U,EJ^ 
coordinate  system  has  its  origin  at  the  DSP 
satellite  making  the  observation  and  the 
directions  of  the  coordinates  are  as  shown. 
Defining 

s  =  {(p,X,hf 

r  =  (x,y,zf  (14) 

\r  =  iU,E,Nf 


T^  =  Tias.HT^-T,as,)/^  (18) 

b)if  7'm.x  <Tnext 

7'6.=7’to  +  (7’„.x-7’««r)/2  (19) 

State  Vector  Generation 

The  state  vector  completely  defines  the 
TBM’s  position  and  velocity.  With  the  tactical 
parameters  and  burnout  time  estimates  the 
altitude  and  range  at  burnout  are  obtained  using 
Eqs.  (1)  and  (2)  with 

tto-T^-To-  (20) 

Using 


A  A 

A  A  -  ^ 

^  dr'^*  dw 


(15)  dbo  =  dbolK 

where  is  the  Earth  radius 


(21) 
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4>ho=  (ft  1 2) -cos 


-1 


'cos0^sin0o  + 
sin  0^^  cos  00  cos  aQ 


^bo  =  ^o+  sin~‘ (sin  9^  sin ag  /cos(l>to)  (22) 

alt  to  =  ho +  hbo 


This  gives  the  position,  the  velocity  at  burnout 
is  expressed  in  terms  of  the  speed  the 

flight  path  angle  (y^)  and  heading 


rto  =  m-^(h/d) 

a  to  =  sin"^  (cos  0o  sin  Kq  /  cos  ^to ) 
d  =  ddldt,k  =  dhldt 


The  state  vector  at  burnout  is 

(^io  Tfco  )  • 


Impact  Position  and  Time 

The  ballistic  trajectory  is  modeled  in 
three  phases  powaed  flight,  which  is  the 
portion  DSP  observes;  free-flight,  which  is  a 
portion  of  an  elliptic  orbit  or  ballistic  trajectory; 
and  re-entry,  which  is  the  portion  in  which 
atmospheric  drag  becomes  significant  until 
missile  impact.  The  powered  flight  is  modded 
with  the  TBM  profile  polynomials  presented 
earlier.  The  ellipse  trac^  during  free-flight  is 
simulated  using  inertial  two  body  mechanics. 
Atmospheric  drag  effects  during  the  re-entry 
phase  are  not  specifically  calculated  in  this 
model,  but  are  somewhat  accounted  for  by 
assuming  that  the  distance  traveled  over  the  Earth 
from  re-entry  to  impact  is  the  same  as  fiom 
launch  to  burnout,  lliis  is  the  same  as  assuming 
that  the  Earth  central  angles  are  the  same,  i.e.. 
These  are  approximations,  but  their 
effect  should  be  negligible  for  the  error  analysis. 
They  would  not  be  negligible  if  one  was 
concerned  with  the  actual  trajectory.  See  Ref.  2 
for  details. 


Error  Sources 

Each  error  source  usually  has  one  or  more 
underlying  causes.  For  example,  the  LOS  aror 
is  a  result  of  attitude  errors  and  focal  plane 
misalignments.  A  complete  error  analysis  would 
break  each  aror  source  down  to  its  fundamental 


level,  and  model  each  level  correctly.  In  this 
section,  some  of  the  underlying  causes  of  the 
overall  sources  are  identified,  but  in  these 
analyses,  only  the  overall  eriOT  magnitudes  ate 
analyzed.  The  errors  are  grouped  into  three  types: 
time,  LOS  and  satellite  position.  Time  errors 
can  be  caused  simply  by  having  more  than  one 
clock  referenced  as  a  source  of  time  measurement 
Imperfect  synchronization  between  clocks’  time 
and  time  passage  rate  are  obvious  otot  sources. 
Time  delays  caused  by  radio  transmission  of  data 
due  to  distance,  atmospheric  refraction,  and 
relative  motion  Doppler  effects  may  add  another 
time  error.  The  magnitude  of  time  errors  is 
relatively  small,  and  getting  smalls  as  time 
measurements  are  being  made  continually. 

LOS  measurement  errors  can  arise  from 
many  sources,  primarily  attitude  uncCTtainties  and 
IR  radiation  measurement  errors.  Any  attitude 
control  system  inaccuracy  effects  are  amplified  by 
the  geosynchronous  altitude.  The  knowledge  of 
the  telescope  alignment  with  the  satellite’s 
reference  fr^e  is  defined  by  the  design  and 
manufacturing  of  the  DSP  satellite,  and  changes 
slightly  with  thermal  variations.  The  IR 
radiation  measurements  of  intensity  and  angle  of 
arrival  (AOA)  also  have  several  underlying  error 
sources.  Locations  of  the  individual 
photoelectric  cells  on  the  focal  plane  array  ate 
recerded  in  what  is  tamed  the  “E’bcal  Plane 
Vector  Table”  (EPVT).  The  positions  of  the 
PECs  change  as  the  satellite  heats  and  cools  with 
varying  sun-satellite  orientations,  causing  a 
warping  of  the  focal  plane,  but  the  FPVT  does 
not  account  for  the  changes  real-time  .  to 
addition,  detector  noise,  refraction,  and  IR 
attenuation  due  to  clouds  and  water  v^x)r  all  add 
to  the  total  measurement  error. 

'  Satellite  position  measurements  fiom  the 
AFSCN  are  u^  to  update  the  ephemeris  once 
per  week.  The  ephemeris  is  propagated  from  the 
time  of  update  to  estimate  the  satellite  position 
for  the  next  week.  The  measurements  are,  of 
course,  inexact,  and  the  resulting  ephemeris  error 
grows  with  time  during  the  one  week  prediction 
interval. 

Since  the  missile  profile,  Eqs.  (1)  and  (2),  is 
an  approximation  of  an  average  trajectory  for 
each  missile  class  it  is  a  source  of  error.  The 
profiles  can  be  a  major  error  source.  HbwevCT, 
this  analysis  does  not  OHisider  this  mor  source 
and  assumes  the  profile  to  be  perfect 

Another  source  of  error  is  the  burnout  time. 
This  error  is  not  the  result  of  an  error  in  the 
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system,  but  is  the  result  of  the  10  second  scan 
period.  The  error  in  the  burnout  time  is 
considered  in  this  analysis. 

Keeping  all  these  underlying  components  in 
mind,  the  errors  are  modeled  in  the  tactical 
paramet^  algorithm.  The  real-world  oror 
component  magnitudes,  bias  and  random 
distributions  may  be  different  from  those 
simulated,  but  the  overall  effects  manifest 
themselves  in  a  manner  close  to  the  error  models. 
It  should  be  possible  to  extrapolate  the  results 
obtained  in  iis  study  to  diJEfaent  errors  by 
properly  scaling  the  different  magnitudes  and 
distributions.  It  must  be  emphasized  that  since 
the  goal  of  this  study  is  not  to  exactly  determine 
the  tactical  parameters  at  launch,  state  vector  at 
burnout,  or  impact  time  and  position.  The 
purpose  is  to  determine  the  contribution  of  the 
various  error  sources  upon  these  values,  and  this 
can  be  done  with  approximate  models  for  the 
errors.  If  the  intended  goal  is  to  calibrate  the 
system  to  eliminate  bias  errors,  then  the  otot 
sources  must  be  modeled  more  exactly  to 
determine  what  is  observable. 

Results 

The  algorithm  previously  described  was 
programmed  into  MATLAB™.  Three  entr 
sources,  time,  satellite  position  and  LOS,  woe 
simulated  and  added  to  the  observational  data, 
first  separately,  and  then  combined.  The  effects 
of  the  errors  on  the  results  are  analyzed  at  three 
points:  launch  position,  burnout  position,  and 
impact  position.  Five  Middle  Eastern  Capital 
cities  were  chosen  as  fictitious  launch  sites  to 
determine  the  effects  of  TBM  launch  position  on 
the  accuracy  of  the  results.  Except  when  noted 
there  wCTe  three  DSP  satellites  in  geostationary 
orbit  at  longitudes  10°E,  70*’E  and  105'’E.  The 
approach  was  Monte  Carlo  with  1000  runs  fw 
each  condition.  .  . 

The  model  for  the  time  error  is  a  random 
uniform  distribution  between  zcto  and  one 
milliseconds.  For  satellite  position,  a  random 
normal  distribution  with  zero  mean  and  standard 
deviation  of  200  meters  was  added  to  eadi 
components,  radial,  in-track  and  cross-track. 
This  is  approximately  a  one  standard  deviation 
sphere  of  346  meters.  The  LOS  error  was 
modeled  as  a  random  normal  distribution  with  a 
standard  deviation  of  five  microradians,  added  to 
both  the  P  and  ti  components.  The  values  used 


for  the  error  sources  are  representative  of  what 
can  occur. 

Figures  7-9  show  for  a  due  east  launch  fiom 
Raghftvl  the  distribution  and  error  ellipse  for  the 
launch  point,  burnout  point  and  impact  point, 
respectively.  The  error  ellipse  (ellipsoid)  is  the 
figure  defined  by® 

z^p-‘z  =  l  (24) 

where  Z  is  the  vector  representing  the  variables 
in  question,  e.g.,  latitude  and  longitude,  and  P  is 
the  covariance  matrix  for  those  variables.  Recall 
that  the  time  error  and  the  burnout  time  are  not 
Gaussian  so  the  covariance  matrix  does  not 
represent  a  Gaussian  oror  distribution.  Figure  8 
is  an  ellipsoid  since  the  burnout  point  is  tiiree 
dimensional.  Shown  in  the  figure  are  the  three 
orthogonal  ellipses  for  each  set  of  two  axes. 

The  data  for  the  300  km  TBM  were  collated 
by  city,  heading  and  error  source.  The  maximum 
area  (volume)  of  the  error  ellipse  (ellipsoid) 
obtained  from  the  1000  runs  for  each  of  the  three 
points  are  shown  in  Figures  10-18  for  each  of  the 
three  error  sources.  Figs.  19-21  show  the 
combined  error  effects  for  the  three  points.  Each 
figure  has  five  curves  for  the  five  launch  points 
(cities).  Immediately  apparent  is  that  the  launch 
point  has  very  little  effect  on  the  size  of  the  oror 
ellipses  and  ellipsoids.  Also,  evident  is  that  the 
observational  geometry  and  TBM  heading  have  a 
large  effect  on  the  ellipsoid  size. 

Table  2  summarizes  the  results  by  showing 
the  minimum,  mean  and  maximum  area 
(volume)  for  each  error  source  for  the  three 
points.  It  is  evident  from  the  results  that  the 
error  sources  can  be  tanked  in  order  of  effect: 

1.  LOS  errors 

2.  Satellite  position  errors 

3.  Time  errors 

Taking  advantage  of  the  ease  of  modifying  the 
MATLAB™  code,  various  changes  to  the  present 
DSP  system  model  woe  made  to  determine  the 
effects  upon  the  accuracy  of  the  results.  For 
these  cases  the  launch  site  was  Baghdad  and  the 
observational  data  were  modified  by  the  combined 
error  sources. 

The  first  case  is  the  control  case  with  nominal 
parameters  and  is  denoted  by  “Bagcom”  to 
represent  “Baghdad  combined  errors”,  and  is  used 
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Table  2  Error  Ellipse  Areas  and  Volumes 


Error 

Impact  Ellipse  Area 
(km*) 

min 

mean 

max 

min 

mean 

max 

min 

mean 

max 

time 

3.26e-8 

2.59e-7 

1.22e-6 

2.9e-ll 

3.6e-10 

1.33e-9 

1.57e-6 

2.61e-5 

1.41e4 

iRjjSSjSIIII 

3.77e-3 

5.27e-3 

8.68e-3 

6.04e-4 

1.14e-3 

2.32e-3 

6.12e-3 

7.92e-3 

327e-l 

LOS 

2.97e-2 

8.45e-2 

1.83e-l 

2.37e-2 

9.33e-2 

2.09e-l 

3.05e0 

1.02el 

1.95el 

comb. 

3.91e-l 

l.Ole-1 

2.97e-l 

4.44e-2 

1.25e-l 

4.50e-l 

3.35e0 

1.38el 

4.06el 

Tables  Various  Case  Comparison 


Emx 

Launch  Ellipse  Area 
(km*) 

Burnout  Ellipsoid  Vol 
(km*) 

min 

mean 

max 

min 

mean 

max 

min 

mean 

max 

Bagcom 

4.21e-2 

l.Ole-l 

2.74e-l 

4.46e-2 

2.60e-l 

3.83e0 

1.40el 

4.04el 

4.40e-2 

7.07e-2 

l.lOe-1 

2.09e-2 

4.89e-2 

8.61e-2 

3.94e0 

7.82e0 

1.52e0 

Molniya 

4.39e-2 

7.58e-2 

1.25e-l 

2.61e-2 

5.57e-2 

1.04e-l 

2.91e0 

7.96e0 

1.85el 

10  s  sc 

4.21e-2 

l.Ole-1 

2.74e-l 

4.46e-2 

1.09e-l 

2.60e-l 

3.83e0 

1.40el 

4.04el 

7.5  s  sc 

2.62e-2 

7.44e-2 

1.95e-l 

2.53e-2 

4.55e-2 

1.09e-l 

4.03e0 

8.72e0 

2.08el 

5.0  s  sc 

2.48e-2 

5.02e-2 

1.52e-l 

1.23e-2 

2.52e-2 

6.52e-2 

2.52e0 

5.84e0 

1.64el 

2.5  s  sc 

1.18e-2 

2.36e-2 

3.48e-2 

5.05e-3 

9.21e-3 

1.31e-2 

1.03e0 

2.95e0 

4.58e0 

as  a  baseline  case  for  comparison.  The  second 
case,  “Synchr.”  shows  the  effect  of  synchronizing 
the  spins  of  the  satellites  so  that  the  satellites 
scan  the  area  of  interest  within  one  second  of 
each  other.  In  the  third  case,  “Molniya”,  a  third 
satellite  at  TO^E  GHA  and  63.4°  latitude  rq)laced 
the  DSP  satellite  at  the  same  longitude  to 
simulate  a  Molniya  +  geostationary  viewing 
geometry.  The  purpose  hwe  was  to  evaluate  the 
effect  of  better  triangulation  of  the  launch.  The 
remaining  cases  simulated  faster  scan  periods 
from  10  seconds  down  to  2.5  seconds.  The 
results  are  shown  in  Table  3. 

The  first  two  modifications  have  the  effect  of 
decreasing  the  areas  and  volume  by  about  30%. 
A  decrease  is  expected  for  the  Molniya  case  since 
two  satellites  on  the  equator  and  one  at  high 
latitude  gives  better  viewing  geometry.  The 
synchronized  spin  results  were  somewhat  of  a 
surprise,  however.  The  effect  on  the  launch 
ellipse  area  is  fine,  but  since  the  burnout  time 
estimate  should  be  less  accurate,  it  would  be 
expected  that  the  burnout  ellipsoid  volume  and 
impact  ellipsoid  area  would  increase.  Howcvct, 
they  decreased.  More  research  is  needed  in  this 
area. 


Increasing  the  scan  rate  has  two  positive 
effects;  more  data  is  obtained  during  the  boost 
phase  allowing  for  a  better  trajectory  estimate, 
and  a  better  estimate  of  the  burnout  time  is 
obtained.  The  decrease  in  the  launch  ellipse  area 
and  the  impact  ellipse  area  are  essentially  linear 
with  decrying  scan  rate  and  the  burnout 
ellipsoid  volume  decrease  is  quadratic. 

Conclusions 

The  TALON  SHIELD/ALERT  state  vector 
estimation  algorithm  was  modeled  in 
MATLAB™.  System  errors  woe  modeled  and 
introduced  into  the  algorithm  to  determine  the 
effects  upon  accuracy  of  the  final  results.  The 
errors  were  broadly  categorized  as  errors  in  time, 
satellite  position  and  LOS.  The  relative 
magnitudes  of  the  error  effects  listed  largest  to 
least  are: 

•  LOS  •  using  a  zero  mean,  S  |jradian  standard 
deviation  normal  errOT  distribution  in  focal 
plane  coordinates 

•  Satellite  position  -using  a  zero  mean,  200 
meter  standard  deviation  normal  oror 
distribution  in  satellite  position  coordinates 
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•  Time  -  using  a  uniform  error  distribution 
between  zero  and  one  millisecond 
As  expected,  the  effects  behave  independently  and 
the  principle  of  superposition  can  be  used. 


Ai-— ,A2-^ 

dvr  ^  _  dy 


dr 


dvr 
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Appendix 

In  this  Appendix  the  partial  derivatives  are 
presented.  From  Eqs.  (4)  and  (14)  -  (16) 
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Figure  6  Howchart  Diagram  of  Iteration  Process 
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Figure  7  Launch  Position  Error  Ellipse 


Figure  8  Burnout  Position  Error  Ellipsoid 
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Figure  9  Impact  Position  Error  Ellipse 
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Figure  10  Time  Error  Effects  Upon  Launch 
Ellipse  Area 
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Figure  11  Time  Error  Effects  Upon  Burnout 
Ellipsoid  Volume 
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Figure  13  Satellite  Position  Error  Effects  Upon 
Launch  Ellipse  Area 


Figure  14  Satellite  Position  Effects  Upon 
Burnout  Ellipsoid  Volume 


Figure  12  Time  Error  Effects  Upon  Impact 
Ellipse  Area 


Figure  15  Satellite  Position  Upon  Impact 
Ellipse  Area 


Heading  (degrees) 


Heading  (degrees) 


Figure  16  LOS  Error  Effects  Upon  Launch 
Ellipse  Area 


Figure  19  Combined  Error  Effects  Upon  Launch 
Ellipse  Area 
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Figure  17  LOS  Error  Effects  Upon  Burnout 
Ellipsoid  Volume 


Figure  20  Combined  Error  Effects  Upon 
Burnout  Ellipsoid  Volume 
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Figure  18  LOS  Error  Effects  Upon  Impact 
Ellipse  Area 


Figure  21  Combined  Ettot  Effects  Upon  Impact 
Ellipse  Area 
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